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Abstract 

This paper presents new measurements in a case study of the minimal SUGRA 
model with m = 100 GeV, roi/j = 300 GeV, A = 0, tan/3 = 2, and sgn/i = + 
based on four-body distributions from three-step decays and on minimum masses in 
such decays. These measurements allow masses of supersymmetric particles to be 
determined without relying on a model. The feasibility of testing slepton univer- 
sality at the ~ 0.1% level at high luminosity is discussed. In addition, the effect of 
enlarging the parameter space of the minimal SUGRA model is discussed. The direct 
production of left handed sleptons and the non-observation of additional structure in 
the dilepton invariant mass distributions is shown to provide additional constraints. 
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1 Introduction 



If SUSY particles exist at the TeV mass scale, they will be produced at the LHC with large rates, 
so discovery of their existence will be straightforward. In minimal SUGRA and similar models, 
however, the decay products of each SUSY particle contain an invisible lightest SUSY particle 
(LSP) Xii so 110 masses can be reconstructed directly. Previous studies of SUGRA models have 
concentrated upon extracting information from kinematic end points measured in three-body final 
states [§, [|, |], ||, H resulting from decays of the type X2 ~^ Xi^ + £~~ an d q_L —> X2Q ~^ hxiq. Studies 
of GMSB |§| models have shown that multi-step decay chains such as X2 ~~ * ^r^ T ~* X?^ + ^~ — * 
G^jt + l~ can provide multiple constraints which can be used to extract masses without fitting to 
any underlying model. 

In this paper we exploit multi-step decays for SUGRA models, specifically the decay 

QL - X°2Q - ^ T - X°i£ + rq (1) 

for the minimal SUGRA model with the previously studied ||, g] parameters mo = 100 GeV, 
m i/2 = 300 GeV, Aq = 300 GeV, tan/3 = 2, and sgn^ = +. The masses for this point are given 
m Table @. As we shall show below, we are able to reconstruct both upper edges for the £ + £ , 
£ + £~q, and £^q mass distributions and a lower edge for the £ + £~q mass coming from backwards 
decays of the X2 i n the q^ rest frame. (The use of analogous upper and lower edges to reconstruct 
masses has been extensively discussed for e + e~ colliders ||.) These measurements make it possible 
to reconstruct all the masses involved in the decay. As part of this analysis we develop a fitting 
procedure to make estimates of the statistical errors of the various measurements that are more 
quantitative than previous ones (2|, ||. 

We then illustrate how these and other techniques can be used to over-constrain the parameters 
of the minimal SUGRA model and place constraints on the model itself. As part of this study 
we will show how some signals change qualitatively as one varies the relations between the squark 
and slepton masses that the minimal SUGRA model predicts. We allow for the masses of the 
third generation squarks and sleptons to vary and for the masses of the sfermions in the 5 and 
10 representations of SU(5) to differ. This is the first step in assessing how well more general 
supersymmetric models can be constrained using data from the LHC. 

The analyses presented here are based large samples of events (the actual numbers are given 
below) generated using ISAJET [pi] and RUNDST 5, which implements a simple detector simu- 
lation 3 including efficiencies representative of the ATLAS detector. Jets were found using a 
fixed cone algorithm of size R = 0.4. Missing energy was determined including the rj coverage and 
a Gaussian approximation to the energy resolution of the calorimeter. Leptons energy resolutions 
were also included (for details see Ref. [||) and an appropriate detector efficiency of 90% was 
included. The event selection cuts make the Standard Model background small compared to the 
SUSY signal, so that clean SUSY samples can be studied. In an actual experiment, the cuts are 
likely to be less severe so that more signal events will survive. The dominant background with 
our cuts is combinatorial background in the interesting SUSY events and other SUSY events that 
happen to pass the cuts. 

This paper contains treats several distinct but closely related topics. Section || describes the 
extraction of combination of masses from four-body kinematic limits. Section ^ describes the use of a 
lower edge to determine another combination of masses, and Section ^| combines these measurements 
to determine several masses of SUSY particles without relying on a model fit. Section || discusses 
the errors that could be achieved on the dilepton endpoint. We already know |^, |(| that this error 
is quite small, but as precise a measurement as possible is useful as a test of e/fi universality. 
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Table 1: Masses of the superpartners, in GeV in the default case (Point 5) and in two modified 
cases with ms = 75 GeV and = 125 GeV described in Section ||. The first and second generation 

of squarks and slep tons are degenerate and so are not listed separately. 

Superpartner default m§ = 75 GeV = 125 GeV 
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Section || explores additional signatures that are relevant to some simple extensions of the minimal 
SUGRA model. Section |7] reexamines the global fits of parameters both for the minimal SUGRA 
model and for its extensions. Finally, after some concluding remarks, the appendix describes an 
unsuccessful attempt at full reconstruction of SUSY events. 

2 Information from four-body decays 

In this and the following two sections (||, ||) and in the Appendix, the analysis is based on a sample 
of 10 6 SUSY events generated with ISAJET 7.32. This sample corresponds to approximately 70 
fb _1 of integrated luminosity. This large sample is needed so that the the statistical fluctuations 
shown on the plots in the these sections corresponds approximately to those expected in the actual 
data from one year at the LHC design luminosity. The Standard Model background is not included 
on the plots (generating comparable statistics for this is a prohibitive task) but is known to be 
small in the channels 0, |6| that are used in these sections. 
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Previous work on measurements in SUGRA models relied mainly on end points measured in 
three-body final states with one invisible particle || ||, |4|, ||, |6j. It was subsequently found || 
in studying GMSB Point Gla that four-body distributions from multi-step decays such as x.2 ~^ 
£%£* -» jftPlr -» G~f£ + £~ contain even more information. This method has considerable gener- 
ality as we now illustrate by applying it to the decay chain given in Equation [l] above. The three 
observed particles, £ + , £~ and q (which appears as a hadronic jet) can be used to make several 
mass distributions: £ + £~q, £q and £ + £~ . The last distribution was considered previously [||, || and 
has a sharp kinematic end point that results when the unobserved Xi momentum is minimized in 
the rest frame of £ R . 

In order to ensure a clean sample of SUSY events the following event selection was applied. 

• At least four jets with pr,\ > 100 GeV and pt,2,3,4 > 50 GeV, where the jets are numbered in 
order of decreasing px- 

• M e g > 400 GeV, where M e g is the scalar sum of the transverse momenta of the four leading 
jets and the missing transverse energy: 

M e fi =PT,1+ PT,2 + PT,3 + VTA + $T ■ 



• $ T > max(100GeV,0.2M eff ). 

• Two isolated leptons of opposite charge with px > 10 GeV, \r]\ < 2.5. Isolation being defined 
so that there is less than 10 GeV of additional transverse energy in an R = 0.2 cone centered 
on the lepton. 

With these cuts the Standard Model background is negligible, as can be seen from Figure 26 of 
Reference ||. Thus, Standard Model backgrounds will not be shown here. 

It is expected that the two hardest jets will be those coming directly from qx — ► X2Q as a 
dominant production procees is that which leads to q~x9 and hence to pairs of q~x- Therefore, the 
smaller of the two masses formed by combining the leptons with one of the two highest pt jets 
should be less than the four-body kinematic end point for squark decay, namely 



M max 



M| - M 2 ( M? - M? 

v 1L x%) \ X\ X\ 

X°2 



1/2 

= 552.4 GeV 



The distribution of the smaller £ + £~q mass is plotted in Figure |T| with same-flavor lepton pairs 
weighted positively and opposite-flavor ones weighted negatively. The e + e~ + ^ + fJ.~ — e^p^ com- 
bination cancels all contributions from two independent decays (assuming e-[i universality) and 
strongly reduces the combinatorial background. This distribution should vanish linearly as the end 
point is approached. Figure [l] also shows a linear fit near the end point. The extrapolation of 
this fit gives an end point of 568.0 GeV, 3.4% above the nominal value. The distribution itself is 
quite linear, so varying the interval over which the fit was made produces only small changes in the 
end-point value. 

An additional selection was then made: one £ + £~q mass was required to less than 600 GeV and 
the other greater, so that the assignment of the jet to combine with the lepton pair is unambiguous. 
The combination with the smaller mass is then used for further analysis. The mass distribution 
of the £^q sub-system was then calculated for each lepton and the selected jet. If the jet mass is 
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Figure 1: Mass distribution for the smaller of the two £ + £ q masses showing a linear fit near the 
four-body end point. 



neglected, then the mass for the jet and the first lepton emitted has an end point analogous to the 
£ + £~ one at 




= 479.3 GeV 



In order to make the former structure as clear as possible, the combination with the larger invariant 
mass out of the two possible £^q pairings in the £ + £~ q combination is used in making Figure ^. 
Again events entered the histogram weighted by flavor (+1 for e + e~ and events and — 1 for 

e ± /x =F ) in order to reduce combinatorial background. 

There is also an end point where the spectrum vanishes for Mi„ formed using the lepton origi- 
nating from the last step in the decay chain, Equation [j], at 




= 274.5 GeV . 



The structure from this end point is buried under former distribution and hence is not visible. (In 
the GMSB case studied previously ||, the analogs of both end points were visible.) 

If the resolution were perfect and there were no selection cuts, the £q mass distribution for the 
lepton emerging from the second step in the decay chain would be given by 

M ^ ax \l^-^dz, -Kz<l 

where M^ ax is the end point given above and z = cos 9* is the cosine of the decay angle of the 
slepton in its rest frame. In order to determine whether the selection cuts or resolution are more 
responsible for the distortion and to estimate how well this end point might be measured, this form 
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Figure 2: Distribution of the larger of the two l^q masses for i + i~ q events in which Mu q < 600 GeV 
and a fit described in the text. 
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Figure 3: Distribution for the smaller of the two £ + £ q masses from Figure [l]; a Gaussian-smeared 
fit plus a linear background described in the text is also shown. 



was smeared with a Gaussian in an attempt to parametrise the resolution effects. The fit function 
is 

r 1 



r+i 

f(M) = J dzAexp 



with parameters A, M?J ax , and a. The integral was done numerically using 96-point Gaussian 
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quadrature integration, and the fit to Figure ^ was made using MINUIT and MINOS flll] ] inside 
PAW. The resulting fit, shown in Figure g, gives M £ ™ ax = 433.2jJ;a GeV, which is 9.6% lower 
than the true position and a = 58.2 GeV. It has a reasonable x 1 , indicating that the cuts do not 
significantly distort the shape of the distribution over the fitted range. The shift to lower values is 
due primarily to energy lost out of the R = 0.4 jet cone. If the analysis is repeated with R = 0.7 
jet cone it is about 3.8% low. The resolution is due mainly to the resolution on the jet energy 
measurement and is consistent with that expected given the form of our detector simulation. 

The resolution smearing will also shift the position of the Uq end point. This distribution was 
re-fit using the empirical form 



f(M) 



M, 



Uq 



dz (ai{M Uq - z) + a 2 {M Uq - z) [ 



exp 



1 

2? 



(M - zf 



+ h + b 2 M 



using the same a obtained above and fitting for a^, bi, and Mu q . This fit is shown in Figure p] and 



gives M Uq = 498.0^4 GeV > which is 9.8% low. The R = 0.7 jet cone gave a value that was 4.7% 
low. In an actual experiment these shifts due to energy loss out of the cone could be understood by 
using detailed comparisons of Monte Carlo simulations with data and using jets of known energy 
to set the jet energy scale Ji~2]j . 

The ratio of the £q and Uq end points is independent of M qL : 



Mj 



iq 



M Pl 



\ 



Mi 



ri 



Mi 



M\ - M% 



0.868. 



compared with a fitted value of 433.2/498.0 = 0.870 for the fits obtained using resolution smearing. 
This ratio should be less sensitive to the jet energy scale and so measured more accurately than 
the individual end points. The difference between the fitted and computed values is very small and 
the result is stable; repeating the same analysis using jets defined with an R = 0.7 cone shifts the 
individual fitted edges by ~ 5% but gives 0.877 for the ratio. 



3 Lower edges 

The upper limit of kinematic distributions has been used in the previous section to extract infor- 
mation. Kinematic distributions can also have lower limits. These have been exploited for example 
in the NLC SUSY analysis |J. For a process like e + e~ — > fi. + fl~ — > Ai + Xi/ i_ X?) the fixed center of 
mass energy and resulting fixed momentum of the \x results in a maximum and minimum energy 
of the observed muon which can be used to extract both fi and x? masses. 

A similar analysis can be used for the process given in Equation [l] at the LHC. The squark mass 
plays a role analogous the center of mass energy in the e + e~ case and a Lorentz invariant quantity 
must be used. For a given value of z = cos 6*, the decay angle of the second lepton in the res t 
frame, the £ + £~ mass is determined to be 

mi = (Msn 2 ^ , 



where 



M$ ax 



\Ml-Ml){Ml-Mp 
Ml 
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There is a corresponding momentum pu in the x% res ^ frame. Thus as a function of z there is a 
minimum of the Mggq mass. For z = the expression for this minimum simplifies to 



( M mirn2 



L Uq 



M mm 
1V1 Uq 



4M'£M2 X 

-MfM| + "iMlMlMl - M$M 2 - M%M^ - M?M$M%- 
MlM 2 e M 2 q + 3M 2 2 M 2 M 2 - M e 4 M ? 2 + (M| - M 2 )x 
^/(M 4 + M e 4 )(M| + M 2 ) 2 + 2M 1 2 M 2 (M 2 4 - 6Mf M e 2 + M e 4 ) 
271.8 GeV 



(2) 



where 



M 2 



M^o, M e = Mi , Mi = M^o 



In order to extract M^ q n , events were selected as before with the additional requirement Mu > 
M^ ax /\/2, corresponding to z > 0, and the larger of the two possible itq masses formed by 
combining the lepton pair with the two highest px jets was chosen. This distribution is shown 
in Figure ^ The lower edge is not very sharp, presumably because gluon radiation can carry off 
energy and so give masses below the nominal end point. Nevertheless, the lower edge is clearly 
visible and is not obscured by the kinematic cuts. 

A fit to the shape including a Gaussian resolution where the width was allowed to float is un- 
stable; as the fit range was changed the value of the width changed significantly and was sometimes 
unreasonably large. This effect seems to be due to the small tail below the edge. Therefore, the 
Gaussian width was constrained to be 10% of the edge value (M]°™). This width was used to smear 
the form 



[A(M - MjT) + B(M - Mfj%Y\6{M - M t 



low\2i 



liq ) 



which was then fitted to the distribution allowing A, B and M l $% to float. The fit gives = 
283.7i£i GeV - The X 2 for the fit is rather poor (30 for 11 degrees of freedom), mainly because of the 
few bins around 200 GeV; a better fit can be obtained by restricting the range to 300 — 600 GeV. 
Because of this, and because the edge is not very sharp, more study with different choices of 
the SUSY parameters is needed to understand the actual error on the fitted value that could be 
achieved. We shall assume conservatively in the discussion below (see Section that an error of 
±2% can be achieved. 

There is also a minimum value of the hq invariant mass from the decay chain — > X2Q ~^ Xi^Q- 
If the jet is again treated as massless, this is given by 



(Mir 



l hq 



1 

2MZ 



(M t 



Mi) 



346.5 GeV 



Ml) 



(Mf 



Ml 



Ml) 2 - AMlM 2 



(3) 



It has been shown previously that the Higgs decay to bb || ||] can be extracted. The following 
cuts were applied 

• M eS > 400 GeV; 

• $ T > max(100GeV,0.2M eff ); 

• at least four jets with px > 50 GeV and one with px,i > 100 GeV; 
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Figure 4: Distribution of the larger of the two llq masses after cuts described in the text, showing 
the lower edge. 



• Transverse sphericity St > 0.2. 

In addition, events were selected to have exactly two b jets with pT t b > 25 GeV, 70 < M&6 < 
110 GeV, and the larger of the two masses formed by combining this pair with either of the two 
hardest jets was selected, since at least one of the two should be greater than the minimum mass. 
This distribution is shown in Figure [|. While there is a threshold in roughly the right place, it is 
not very distinct. This probably is due to a combination of resolution for this multijet system and 
of the substantial combinatorial background under the Higgs peak. A side-band subtraction and 
careful jet energy calibration might be able to clean up the distribution. We will assume below an 
error of ±5%; this is substantially larger than the llq error and will only add a very weak additional 
constraint. It would be important to study this further in cases where the decay to sleptons is not 
available such as "Point 1" and "Point 2" of Refs. | and |. 



4 Model-independent masses 

In this section we discuss how the measurements in the previous section and those discussed pre- 
viously [§, |(| can be used to determine the masses of the SUSY particles without reference to 
the underlying SUGRA model. The identification of the decay chains is needed but these are 
based on much weaker assumptions. We have a number of measurements all related to the process 
QL - - Ir - X°v 



M « ax 



(Ml-M\){M^-M\) 



1/2 



108.9 ± 0.11 GeV (see Ref §) 



ji/rmax 
• 1V1 Uq 



(Ml L -Ml )(Ml -Ml ) 



M 2 



1/2 



552.4 ±5.5 GeV 
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Figure 5: Larger of the two bbq jet masses, showing the hq mass threshold. 



^ 



R 



M 



1/2 



479.3 ± 5.5 GeV 



Mft (Equation |) = 271.8 ± 5.4 GeV 



There are also two measurements related to x% ~~ * 
• The maximum hq mass (see Reference ^ 



(M, 



max\2 
hq ) 



M 2 h + 



Ml - M| g 



M^ 



Ml 



M lv 



4M2M| ; 



2M 2 



which has the value M^ ax 



522.6 ± 5.2 GeV. 



• The minimum hq mass MJ^ n (Equation ||) = 346.5 ± 17.3 GeV 

A large error is assigned to the hq lower edge because it is not sharp and there is a lot of background 
(see above). 

If the lower edges discussed in the previous section are ignored, there are four measurements 

and four unknown masses masses M q ~ L , M^ 2 , Mf and M^o. Nevertheless, for the errors assumed 

there is a one-parameter family of solutions labeled by M^o, with small uncertainties in the other 

xi 

masses for a fixed value of M^o . This remains true even if the errors are substantially reduced. 

If the lower edges are included, then all four masses can be determined. The errors were 
estimated numerically as follows. The q~L, X2> an d @R masses were generated uniformly within ±50% 
of their nominal values, and the Xi mass was calculated using the dilepton edge (M«), which has 
a much smaller error than the other measurements. The x 2 f° r the remaining measurements was 
calculated, and the point was assigned a probability of exp(— x 2 /2). The resulting distributions 
are shown in Figure |(| The resulting errors range from ±12% for the mass of Xi t° ±3% for the 
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Figure 6: Distribution of the Xii £r, x[>, an d <1l masses satisfying all constraints discussed in the 
text. The fitted widths are about ±12%, ±9%, ±6%, and ±3% respectively. 



mass of qi. If the error on Mf^™ were reduced to ±1%, as might be possible with a more careful 
understanding of the systematics, the error on M^o would be reduced to ±7.5%. The errors are 



highly correlated as can be seen from Figures 



8] which show the scatter plots of M^o vs. M^o 



and Mji R vs. M^o. Of course, the errors on the masses are much poorer than those that arise from 



a fit within the SUGRA model (see Section and Ref. 
assumptions. 



]), but they do not involve any model 



5 Dilepton measurement errors 

In order to provide significant constraints on model parameters and tests of the underlying model a 
number of measurements with comparable errors are needed. In general, measuring one combination 
of SUSY masses very precisely is not particularly useful if the other combinations involve jets and 
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Figure 7: Scatter plot of M^o vs. M^o solutions satisfying all constraints discussed in the text. 
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Figure 8: Scatter plot of vs. M^o solutions satisfying all constraints discussed in the text. 

so are only measured with an accuracy of several percent. An important exception is the decay 
X® —> — > Xi^ + £~ ) which has an end point at 



M~ 



M~ 



(4) 



13 



400 




50 100 
M,. (GeV) 



150 



Figure 9: £ + £ mass distribution showing the x 2 MINUIT fit using PAW. 



A difference in the end points for e + e~ and n + fJ-~ would directly indicate a difference in the e~R and 
fin masses, which obviously an important issue for testing models that purport to understand flavor 
physics. This decay is generally allowed in SUGRA models which give cosmologically interesting 
cold dark matter fll3| , such as the one discussed here. It is also common in GMSB models since 
the £r has only U(l) couplings and tends to be light. The derivative of the end-point with respect 
to Mg vanishes at the geometric mean of the Xi an d X2 masses but in general is of order one; for 
the masses in the case studied here, 

A M max 

" 0.478. 



dM iR 



The same sample of 10 6 SUSY events was used to estimate how well such an edge might be 
measured with full LHC luminosity. In the absence of cuts, the mass distribution should be given 
by the same formula as discussed in Section B namely 



(M max )2 i + f 



dz 



with z uniformly distributed. This form was smeared with a Gaussian using numerical integration 
as in Section |2[ Figures ||] and |l^ show the resulting fits using MINUIT with either the x 2 or the 
maximum likelihood method; the parameters are the overall normalization, the end point, and the 
Gaussian width. The fitted end points with errors from MINOS Q are 108.71+°-°^ GeV and 



108.60 



+0.065 
-0.060 



GeV respectively. The fits are consistent with each other, 



0.088 

but neither quite agrees 



within errors with the expected end point at 108.92 GeV. The statistical errors are slightly better 



than the systematic errors expected from the lepton energy scale [12| and are comparable to the 
errors on the W mass expected to be achieved ultimately at the Tevatron [Ojand LEP [|IT|. 

The maximum \r]g\ for the either of the two leptons is plotted in Figure |ll| and peaks around 
?7max = 1- Thus both the barrel and the endcap regions of the detector are important. If precise 
electron and muon measurements were available only in the barrel, about half the events would be 
lost. 
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Figure 10: £ + £ mass distribution showing the maximum likelihood MINUIT fit using PAW. 
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Figure 11: Maximum \rji\ for dilepton events. 



Clearly a lot more work is needed to understand how to calibrate the detector and to extract 
information at this level of accuracy. In particular, the discrepancy between the fitted and calculated 
end points even in this highly idealized simulation needs to be understood, presumably by studying 
many different samples of SUSY events. It seems clear, however, that statistical errors below about 
0.1% may be achievable, especially if the masses were a bit lower so that the X\h decay were 

absent. 
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6 Non-minimal SUGRA models 



The case studied so far assumes that the scalar masses are all equal at the GUT scale, an assumption 
that is rather restrictive and may not be valid. By studying variations in this assumption, we can 
try to estimate howl the various LHC signals are modified and how well this assumption could 
be tested. We shall show that qualitatively new signals emerge in our preliminary study of non- 
universal SUGRA (NUSUGRA) models which is carried out using two kinds of NUSUGRA models 
closely related to the SUGRA case discussed above. 

6.1 Variations of masses with SU(5) representations 

We vary the scalar masses at the unification scale by assuming that squarks and sleptons which are 
in the 10 of SU(5) have a common scalar mass tuiq while those that are in the 5 of SU(5) have a 
mass 7715. Here, 777,10 is kept at its nominal value of 100 GeV while 7775 is shifted. Two points have 
been studied in detail, namely 7775 = 75 GeV and 7775 = 125 GeV. The masses of the superpartners 
at these points are given in Table g. Squark masses are almost insensitive to these changes in 7715; 
the shifts are much smaller than the errors that were obtained in Section |4] because m^ 2 plays a 
dominant role in these masses via the strong coupling of squarks. The significant changes take place 
in the slepton spectrum, in particular in the masses of £l- Samples of 200000 events were simulated 
for each of the new cases; 10 6 Standard Model background events were also used to ensure that the 
cuts are effective in disposing of it. 

6.1.1 The m 5 < 100 GeV case 

In the benchmark case, 771,5 = mio = 100 GeV, the decay sequence x\ ~* ~* Xi^ + £~ is allowed, 
giving a very clear signature: the lepton pair invariant mass distribution presents a sharp edge near 
the kinematic limit, Equation ||. For modified points with 777-5 < 85 GeV, the left handed slepton 
becomes lighter than x® thus the decay sequence X2 ~^ ~^ is also allowed, giving rise 

to a second edge at 



The position of the edge is sensitive to M; and thus to 7775. As there is less available phase space 
in this decay than in the decay to £r, these leptons are softer and it is necessary to lower the cuts 
as much as possible to ensure good acceptance. ATLAS expects to be able to detect muons down 
to pt = 5 GeV, so the following selection cuts are applied: 

• M eS > 800 GeV; 

• E T > 0.2M efr ; 

• at least one R = 0.4 jet with px > 100 GeV; 

• £ + £~ pair with \rji\ < 2.5, px,e > 10 GeV, and px^ > 5 GeV; 

• £ isolation cut: Er < 10 GeV in R = 0.2 around the leptons; 

• Transverse sphericity St > 0.2. 



Mi 



r/max 

U 




34.9 GeV for m 5 = 75 GeV . 
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Figure 12: M« distribution after the cuts at the ms = m\Q = 100 GeV point (dashed line) and 
modified point with m§ = 75 GeV (solid line). 
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Figure 13: M„„ distribution after the cuts at the 7715 
modified point with 7715 = 75 GeV (solid line). 



mio 



100 GeV point (dashed line) and 



Figures 12 and 14 show the lepton pair invariant mass for the benchmark case and the modified 
cases with 7715 = 75 GeV and ms = 50 GeV. Figures [l^ and 15 show the muon pair invariant mass 
for the same cases. The edge at low Mm for 777,5 = 75 GeV is clearer in the muon case due to the 
increased acceptance at low pt- The presence of two structures with comparable rates enables one 
to deduce the presence of two decay chains and to measure the two end points. Notice that as m§ 



17 



200 



2 150 
> 

CD 

o 

1 100 

<x> 
> 

LU 



50 



~i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 



\\\ 



J I I I I I I I I I L 




50 100 150 

M,. (GeV) 



200 



Figure 14: M« distribution after the cuts at the ms = m\Q = 100 GeV point (dashed line) and 
modified point with 7715 = 50 GeV (solid line). 
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Figure 15: distribution after the cuts at the 7715 = mio = 100 GeV point (dashed line) and 

modified point with 7715 = 50 GeV (solid line). 

is reduced to 50 GeV, the higher mass structure is becoming weaker because the y^Uh coupling is 
larger than the %2^R one - ^ s m 5 increases, Mj, increases, and the branching ratio for X2 &l 
vanishes quickly around 1715 = 80 GeV. Hence, = 75 GeV is about the upper limit where one 
can distinguish from the benchmark case using this channel. 
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6.1.2 The m 5 > 100 GeV case 

In this case there is no visible effect on the £ + £~ distribution or any of the other distributions 
studied in the previous sections. Production of £l via the decay of strongly interacting sparticles is 
small, so one must rely on direct production via the Drell-Yan process. It is possible || to extract 
a signal for Drell-Yan production with li — ► j^i by requiring two isolated leptons, missing energy, 
and no jets, but the slepton mass can only be inferred from the rate and kinematic distributions. 
The jet veto is essential to elliminate events with leptons arising from the decays of squarks and 
gluinos. 

If 7715 is somewhat larger than 100 GeV, the £l becomes heavy enough that the decay £l — ► x\t 
is allowed. Then the decay chain 

h + h 
I I 

i 

I 

results in a final state with four isolated leptons. The signature is two same flavor opposite charge 
(SFOC) lepton pairs and no jet activity. In order to select these events, the following cuts were 
applied: 

• no jet with pr > 40 GeV and \t]\ < 5; 

• at least 4 leptons with px > 10 GeV and \rj\ < 2.5 forming 2 SFOC pairs; 

• the invariant mass of at least one of the pairs is less than 109 GeV (so that it is a candidate 
to arise from the decay of xB) 5 

• £ isolation cut: Et < 10 GeV in R = 0.2 around the lepton. 

The invariant mass of the three leptons coming from the same left-handed slepton provides 
information about its mass. The three leptons were selected as follows: 

• a pair with invariant mass smaller than 109 GeV is assumed to come from the X2 decay; 

• the remaining lepton with lowest px is assumed to come from the same left handed slepton 
as the pair. 

The invariant mass of the trilepton system is then computed. As the production rate is small, high 
luminosity is needed and an integrated luminosity of 300 fb^ 1 , representing the ultimate that can 
be achieved at LHC, was assumed. The Mm invariant mass should have an upper limit given by: 



\ 




which is approximately 128.0 GeV for 7715 = 125 GeV. 



Figures 16, 17 and [18| show the tri-lepton invariant mass respectively at point 5, and the cases 
with ms = 115 GeV, and 7715 = 125 GeV. As M% increases, the £l X2^ branching ratio increases 
Vo for point 5, 10% for m$ = 125 GeV). A clear signal appears for 7715 = 115 GeV. We estimate 
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Figure 16: distribution after the cuts for the 777,5 = mio = 100 GeV point. Note that the 

number of events is very small. 
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Figure 17: Mm distribution after the cuts for the modified point with 7775 = 115 GeV. 



that for masses larger than this, the position of the edge can be measured with a precision of 3 
GeV and is very sensitive to 771,5. As 777.5 increases further, the production rate falls off and the 
signal disappears for 777,5 above 250 GeV. 



6.2 Variations of masses for the third generation 

Here we investigate the possibility that ttt-o for the third generation squarks and sleptons is different 
from that for the first two generations. As in the 777.5 cases, the greatest sensitivity to this change 
is in the slepton sector, so we are forced to consider the detection of final states containing taus. 
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Figure 18: M(g£ distribution after the cuts for the modified point with = 125 GeV. 



In another study [16], the use of hadronic tau decays was illustrated. It was shown how the tt 
invariant mass distribution could be inferred from the observed decay products and the kinematic 
end point extracted. This method enables m T ~ R to be constrained. If this method could be exploited 
in the case of interest here, slepton universality could be tested with great accuracy. The method 
fails for several reasons all of which are related to the observable event rate and background. 



£r£, £ being e or /u, are 



The signal is less clear because the channels Xt, 
open here, whereas the x\ decays only to staus for the case studied in Ref. 16. The decay to h 



h Xi and 



v° 

A 2 



also generates tt final states at a comparable rate and distorts the shape of the distribution. 

Because the gluinos and squarks are heavier, the total SUSY cross section for the cases 
discussed here is smaller than in the case studied in Ref [16], resulting in smaller event 
samples for the same integrated luminosity. 



Figure 19 shows the tt invariant mass reconstructed from the visible decay products for our bench- 
mark case after the following selection cuts are applied: 

• $ T > 0.2M cS or > 100 GeV; 

• at least one R = 0.4 jet with pt > 100 GeV; 

• at least three additional R = 0.4 jets with px > 50 GeV; 

• two jets identified as hadronic tau decays according to the methods described in Ref. |l~6| . 

The distribution is shown for the subtracted combination t + t~ — t^t^ as this is reduces the 
background from jets that are misidentified as taus. 

A sample of 100000 events have been generated; the plot is normalized to 10fb _1 , which would 
correspond to corresponding to ~ 230000 events so the statistical fluctuations are somewhat bigger 
than they would be in the actual experiment. The signal, of about 40 events in the bin at 50 
GeV, is completely buried by the statistical fluctuations. Even several years of data taking at low 



21 



50 

H — 

o 

> 
o 
O 
o 

™ o 

03 

s= 

CD 

i5 



-50 - 



~i — i — i — i — I — i — i — i — r 



j I I I i I I I I I L 



100 200 
M TT (GeV) 



300 



Figure 19: M T + T ~ distribution at point 5 after the cuts and after subtraction of same sign pairs. 
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Figure 20: M T + T - distribution at modified point with tanfi = 5 (other parameters at their nominal 
point 5 value) after the cuts and after subtraction of same sign pairs. 



luminosity would not be enough to reduce them to a satisfying level. LHC high luminosity would 
give enough statistics, but the M TT reconstruction method used here is based on full simulation [17] 
and has not been proven to be viable at high luminosity; additional pile-up events could compromise 
it. 

To emphasize the deterioration of the signal due to the x% ~^ h°Xi channel, a model with tan [3 
shifted from 2.1 to 5 (the other SUGRA parameters remaining at their nominal values) has been 
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Table 2: Masses of the relevant superpartners, in GeV, in the default case and at modified points 
with m 3 rd = 30, 70, 150 and 200 GeV. The first and second generation of squarks and sleptons are 
degenerate and so are not listed separately. 

Sparticle default m 3 rd = 30 GeV m 3 rd = 70 GeV m 3 rd = 150 GeV m 3 rd = 200 GeV 
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studied. Since is now heavier, this decay channel is closed and the rate of tau pair production 



is 3 times greater. Figure |20| shows that the excess of t + t~ pairs now becomes visible; there are 
approximately 120 signal events in the plot concentrated in the bin at 50 GeV. 

An alternative method of extracting evidence for excess r-pair production based on leptonic 



final states is now illustrated. The method is similar to that of Ref [18]. In order to illustrate 
its sensitivity, we have studied models in which the SUGRA parameters remain at their nominal 
values except that the third generation squark and slepton masses (namely t[, b r , t r , L\, L r ) are set 
equal to m 3 rd at GUT scale. The masses of the relevant superpartners for several values of m 3 rd 
are given in Table ^, from which it can be seen that the largest effect is in the stau masses: as m 3 rd 
increases the taus masses rise and the branching ratio for ~^ 1S reduced. The channel closes 
for m 3 rd > 200 GeV. 

Samples of 200000 SUSY events were generated in each case shown in the table, (except for 
m 3 rd = 30 GeV and 200 GeV where 100000 events were generated). The Standard Model back- 
ground has been added using a 1 million event sample. The only significant background after 
applying our selections is from tt events. In order to select SUSY events with slepton pair, and 
reject Standard Model background, the following cuts have been applied: 

• M eS > 500 GeV; 

• $ T > max(0.2M eff ,250GeV); 

• at least one R = 0.4 jet with px > 100 GeV; 

• at least four R = 0.4 jets with pt > 50 GeV; 

• 11 pair with px.e > 10 GeV, |%| < 2.5, £ being electron or muon; 
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Table 3: Ratios of production rates for lepton pairs for the five models studied (see text) 



m 3 rd 30 GeV 


70 GeV 


100 GeV (point 5) 


150 GeV 


200 GeV 


roc = e+e e ~ ± + ; T V ~ 2.61 


3.86 


3.99 


4.62 


4.38 


rsc = e±e± e i^ 0.88 


0.79 


1.25 


1.05 


1.01 



• £ isolation cut: E T < 10 GeV in R = 0.2; 

After these cuts the signal exceeds the Standard Model background by more than a factor of 
ten; the background is mainly due to top quark pair production. Lepton pairs in the signal events 
arise mainly from two processes. 

• The decay of X21 y ia X2 ~^ t^P* — ► Xi^ + ^~ where £ is e#, ftR, f l5 or f 2 ; I is the corresponding 
lepton. This channel produces exclusively oppositely (OC) pairs, of the same flavor (SFOC) 
in case of selectrons or smuons, and both same (SFOC) and opposite flavor (OFOC) pairs in 
case of leptons from resulting from leptonic tau decays. 

• The decay of a xt pair, produced in the decays of gluinos: even if the g — > xtlQ branching 
ratio is small, this channel is important since it can produce same charge (SC) pairs, as 
gluinos are Majorana fermions. The flavors of the leptons are uncorrelated. The main xt 
decay modes are: xt ~> XiW^ — * Xi^^i^t) Xi ~* Xi^^e xt ~ > ^ being e, fi or r. 

A violation of e, fj,, r universality will be revealed by comparing the number of events containing 
same flavor lepton pair (SF), with the number of events containing opposite flavor lepton pair (OF). 
Define 

e + e~ + n + n~ 



This ratio decreases as the branching ratio for X2 ~* ^ T increases. The second class of processes 
listed above is sensitive to violations of e\x universality only and consequently rsc defined for same 
charge lepton pairs is independent of violations of e, r universality. 

Tablej3] and Figure |2l| show roc an d rsc f° r different values of m 3 rd The error bars shown on 
Figure |2l|correspond to 10 fb _1 of integrated luminosity, rsc is insensitive to violations of e/fj, 
universality and should be the same for all the cases considered. The errors indicated on Figure |2l| 
show that the apparent differences are statistical fluctuations. 

In the region between the benchmark case and m 3 rd = 70 GeV, the contribution to roc from the 
second class of processes decreases dramatically since, at m 3 ,d = 70 GeV, v T becomes significantly 
lighter than xt^ an d the decay mode xt ~* ^ t t becomes important, with a branching ratio of 30% 
instead of 0.4%, to the detriment of decays to electrons and muons. For m 3 rd ~ 75 GeV, T2 becomes 
lighter than X2 ( see table |2|), and the decay channel X2 ~^ ^2 T , opens resulting in a large violation 
of e/r universality and increasing the number of OFOC pairs. For 75 GeV < m 3 , d < 100 GeV, very 
slight change in the ratios is to be expected, since these channels are not open, and the branching 
ratios of the usual channel X2 ~^ ^ changes only slowly. For m 3 rd > 100 GeV the only interesting 
decay mode is X2 ~^ as increases and tends to rn^o, the decay branching ratio x\ ~^ tt 
closes, and roc increases. The channel closes at m 3 rd = 200 GeV. 
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Figure 21: roc an d r sc f° r various m 3 rd values (see text). The error bars correspond to an 
integrated luminosity of 10 fb _1 . 



To summarize: for m 3 rd < 100 GeV, roc is very sensitive to m 3 rd. If 10 GeV < m 3 rd < 75 GeV, 
one should be able to constrain it with an accuracy of a few GeV. For m 3 rd > 100 GeV, things are 
not so easy: for (m 3 rd < 200 GeV), roc increases slightly, giving some sensitivity to its value. As 
a result, models with m 3 rd above 200 GeV can hardly be distinguished from each other. 

7 Determination of parameters 

In previous studies Q, || we needed to rely on a global fit to a specific model, e.g., minimal SUGRA, 
to determine individual masses from the combinations of masses measured by various end points. 
The analysis in Section [i] allows us to extract masses in a rather general way, but the global fit is 
still useful. Since the measurements in Sections |2| and ^ provide new information, we reevaluate 
here the precision with which the model parameters can be determined. The strategy is the same 
as before: the parameter space of the SUGRA model is searched using random sampling in order to 
determine the ±34% confidence limits resulting from the assumed "experimental" quantities and 
their estimated errors. 

The following quantities and estimated errors are used in the fit: 



(M,7 X ) 2 



Ml 



R 



Ml + (M| - M| 
(552.6 ±40 GeV) 2 ; 



Mi 



\ 



M% 
M% 

*2 



M\ 



Ml 



M k + 



(Mf 

*2 



Ml 



M\ 



2 \2 



4MM 



2Mf 

X° 2 



Mf 



108.9 ±0.1 GeV: 



0.865 ±0.02: 
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(M, 



min N i2 



(M™ 1 ) 2 



\ 



(K- M l°)( M h 2 - M V 



478.1 ± 40 GeV 



4MfM 2 



-MfM 4 + 3M?M$M% - M$M% - M^M* - M?M$M* 



MlM 2 e M 2 + 3Aff Aff Af? - M e 4 M 2 + (Af| - M 2 ) x 



(M x 4 + M e 4 )(M| + M e 2 ) 2 + 2M 1 2 M|(M 2 4 - §M%M% + M e 4 ) 



(271.8 ±5.4 GeV) 2 ; 



2M| 



-(M 2 -M 2 )x 



[Ml + M 2 - M 2 ) - 
(346.5 ±17 GeV) 2 . 



(M| — M 2 — M 2 ) 2 - 4M X 2 M 2 



In addition we include in the fit with an error of ±3 GeV. The experimental error on the mass 
from h — > 77 will be considerably less than this. The error reflects our estimate of the theoretical 
uncertainty in relating the Higgs mass to the parameters of the SUGRA model. 

A fit of the minimal SUGRA model to these inputs results in the following values of the pa- 
rameters: 

• m = 100.0 ± 3.63 GeV, 

• m 1/2 = 300.0 ± 4.99 GeV, 

• tan/3 = 2.11 ±0.18, 

• M = +i; 

The errors are symmetric, unlike the earlier fits. Recall that we previously [§] quoted:- 

• m = 100±| 2 GeV, 



m 



1/2 



300+5 GeV, 



+0.3 
-0.5' 



• tan/3 = 1.8 

• H = +1; 

Thus the new measurements improve the fit to the minimal SUGRA model as well as allowing 
masses to be extracted without assuming the model. 

A completely general model at the GUT scale would have as many parameters as the MSSM. 
To keep the problem tractable, we consider three variants of the SUGRA model, each with only 
one additional parameter. We use the same fitting procedure to estimate how well these additional 
parameters could be constrained if the actual data corresponded to the benchmark case, i.e. we 
estimate how well we can actually constrain to SUGRA model. 

We first allow the values of mo at GUT scale to be different for the two Higgs representations. 



Restricting m H 



m 



m h-GUT > 0, leads to a 5-parameter fit and the following result (mo is 



now the common mass for all the other scalars) 
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• mo = 100 ± 3.68 GeV, 



• m 1/2 = 301 ± 5.94 GeV, 

• tan/? = 2.11 ± .18, 

• m = +i; 

• mh-GUT < 430 GeV, 95% confidence 

The insensitivity to mh-GUT arises because the derived value of /U is large (~ 500 GeV) and the 
value of mHjj (the Higgs doublet that couples to charge 2/3 quarks) at the weak scale is determined 
mainly by the top quark Yukawa coupling and stop mass and not by the value of mn v at the GUT 
scale (mh-GUT)- Therefore the masses of Xv X% an d h are insensitive to mh-GUT unless it is 
very large. If the H and A Higgs bosons (and the heaviest gauginos) could be observed and their 
masses measured mij v could be constrained. The masses of these particles vary by ~ 40 GeV for 
parameters in the allowed range. The production rates for these particles at LHC are extremely 
small and their discovery is probably not possible there. This insensitivity to mh-GUT is quite 
general [|l9| . 

We next split the squark and slepton masses at the unification scale: the particles that are in 
the 10 of SU(5) are assumed to have common scalar mass m\o and those that lie in 5 of SU(5) are 
assumed to have common scalar mass m§. Fitting for these, we get 

• mio = 100 ± 3.8 GeV, 

• mi/2 = 300±7° GeV, 

• tan/? = 2.11 ± .23, 

• m = +i; 

• m 5 < 420 GeV, 95% confidence 

The lack of precise constraints on these new parameters can be understood. Since mi/ 2 is sig- 
nificantly larger than mo, the values of the squark masses at low energy are controlled by m\ii- 
The excellent constraint on m\§ arises because it controls mg R which is very precisely determined 
by the t + t~ end point. Il would be observable if the decay x% ~^ ^L^ - ► £ + £~Xi were open as 



discussed in Section 6.1.1. The failure to observe this decay then constrains me L and consequently 
m$. Adding this constraint excludes the region 7715 < 75GeV. Direct production of right handed 
sleptons excludes the region 7715 > 115 GeV although this difficult search requires high luminosity 



as explained in Section 6.1.2. Using these constraints we infer: 



• m 5 = lOOi^Q GeV. 

Finally, we consider the case where the third generation squark and slepton masses at the GUT 
scale are allowed to vary. As can be seen from Table | , the sensitivity is confined to the stau masses 
so that a fit without the information from Section "x2 results in almost no constraint on m^d- The 
errors on the other parameters are as above. Adding the constraints from this section implies that 



m 3rd = 1007* GeV 
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8 Conclusions 



In this paper we have demonstrated a number of new techniques that could be used to determine 
masses and decay properties of super symmetric particles at the LHC. We have shown that cascade 
decays with several steps can be used to reconstruct the masses of super symmetric partilces with- 
out any knowledge of the underlying model. We have illustrated new signals that appear when 
the SUGRA model is extended to have more parameters and have shown in particular how e/r 
universality could be tested. We have further demonstrated the very high precision with which 
many of these parameters can be constrained. 
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Appendix: Attempt at complete reconstruction of SUSY events 

The three-step decay chain — ► \2 ~ y ~~ * Xi provides three mass constraints using the values 
determined in Section |4| and so a OC fit for the Xi momentum is possible. If the same decay is 
selected on both sides of the event, then, in principle, one could completely reconstruct the event 
using I^t to select the best solution. This method was successful for GMSB models with decays 
involving leptons and photons, but it fails in this case due to the experimental resolution, as we 
will now show. 

Events were selected to be consistent with two X2 ~^ -t ^ T — > Xi^ +j ^~ decays and so to have 
four leptons, at least two jets, and missing energy: 

• M eS > 400 GeV. 

• $ T > max(0.2M eff , 100 GeV). 

• At least two jets with pr t i > 100 GeV, px,2 > 75 GeV and at least two charged tracks in 
# = 0.4. 

• Four isolated leptons with pj- > 10 GeV, rj < 2.5. Isolation being defined so that there is less 
than 10 GeV of additional transverse energy in a cone R = 0.2 around the lepton direction. 

The cut on the charged multiplicity of the jets was made to eliminate electrons and hadronic tau 
decays from the jet sample. It was further required that there be one and only one way to form two 
opposite-sign, same-flavor pairs of the four leptons with 20 GeV < Mu < 115 GeV. This ensures 
an unambiguous paring of the leptons consistent with two X2 decays. 

The events were then fit to the hypothesis that the four leptons and the two highest px jets 
came from the q~L — ► X2 ~ * & ~ * Xi decay chain. For each such decay chain there are three mass 
constraints, namely 

(Px°+Pl,i) 2 = M l 
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Figure 22: difference Ap/p between the generated and the best reconstructed \p\. 
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There are also two constraints from Since the measurement errors on the jets are comparable 
to those on the jet energies were smeared by factors A, distributed in Gaussian manner, and 
the best solution was taken to be the one which minimizes 



X 



Plx - P2x) 2 {$y ~ Ply ~ P2yY 



+ 



A? 



+ 



o\\x) a 2 (A 2 ) ' 



where both the missing energy resolutions a($ x ^y) and the jet energy scale resolutions <r(Ai.2) are 
determined using the Gaussian calorimeter resolution. The resulting difference Ap/p between the 
generated and the best reconstructed \p\ for the Xi is shown in Figure 22. A similar reconstruction 
was successful in the case of the GMSB study [0] where the decay chain 



X°2 



was used. The reconstruction works much more poorly than in the GMSB case; while a few events 
are correctly reconstructed, most are not. This is not very surprising; the GMSB case relied on 
leptons and photons which have much better energy resolution than the jets used here. In addition, 
in this case, the x\ momenta are significantly smaller than the jet momenta and so the errors on 
the jet energy measurements are very important. 
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